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ABSTRACT 

We present a large q expansion of the 2d (/-states Potts model free energies 
up to order 9 in 1/y/q. Its analysis leads us to an ansatz which, in the first- 
order region, incorporates properties inferred from the known critical regime at 
g = 4, and predicts, for q > 4, the n th energy cumulant scales as the power 
(3n/2 — 2) of the correlation length. The parameter-free energy distributions 
reproduce accurately, without reference to any interface effect, the numerical data 
obtained in a simulation for q = 10 with lattices of linear dimensions up to L = 50. 
The pure phase specific heats are predicted to be much larger, at q < 10, than the 
values extracted from current finite size scaling analysis of extrema. Implications 
for safe numerical determinations of interface tensions are discussed. 
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Much effort has been recently devoted to the 2-d Potts model, both with 
numerical and analytical techniques. In the former case, the goal was either to 
test numerical algorithms and criteria for distinguishing first-order from continuous 
transitions, or to learn how to extract previously unknown quantities such as 
the interface tension. Although much progress was accomplished, there remains 
some unsatisfactory issues such as, for example, slight inconsistencies in finite 
size scaling analysis of the energy cumulants close to the transition temperature 
/3 t _1 , and discrepancies between exact results and numerical simulations for the 
interface tension. This question is important since only numerical simulations can 
determine this quantity in other cases of physical interest such as the 3-d q = 3 
Potts model or QCD at the deconfinement transition. 

Analytical works have shown [1] that close to flt, the partition function Z of 
the Potts model, in a box of volume V = L 2 with periodic boundary conditions 
(used all through in this work), is equal to the sum of the 'partition functions' 
Zi of the q + 1 pure phases, and a term that falls of exponentially faster in the 
linear size of the system. We shall presently ignore this latter term which contains, 
amongst others, the interface tension effects, and concentrate on the i th phase free 
energy Fj = InZi/V, which is V independent and differentiable many times with 
respect to the inverse temperature j3 at j3 t . 

In this letter, we construct explicit formulae for the ordered and disordered 
free energies F Q and F d of the 2-d (/-states Potts model at q > 4. At large q we 
perform their expansion in power of 1/y/q. At low (q — 4) we conjecture their 
behaviour from their known critical behaviour at q = 4. We show that these two 
descriptions match in a large intermediate q region. Then, using the above rigorous 
result, we add up the Z^s so obtained and make absolute predictions on energy 
probability densities in very good agreement with numerical data. We show that 
the difficulties encountered in the finite size scaling analysis of numerical data are 
due to very large high order cumulants, which has consequences for the extraction 
of interface tensions. 

Many properties of the model are known exactly [2]. In particular, it exhibits 
a temperature driven phase transition which occurs at (3 t = ln(^/g + 1). The 
transition is second-order for q < 4, and its critical properties are described, e.g., 
by the a and v indices, which at q = 4 take the common value 2/3. Accordingly, 
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the correlation length and the specific heat there diverge as 



&=4~i/?-a r 2/3 , c Q=4 ~i /3- a r 2/3 . (i) 

Hence in the vicinity of (3 = /?t, the ratio C/£ remains finite at q = 4_. 

The first-order transition region is q > 4. There the energies E Q and -Ed of 
the ordered and disordered phases respectively are exactly known at (3 t [3] . Recent 
works [4,5] on the largest correlation length at flt have shown a common behaviour 
as q — > 4_|_ 

1 7T 2 

f= — ^x(l + C(x" 2 )) with x = exp( =- - ). (2) 

Sv^ 1 " Pl 21ni( v ^+ v ^4) J lJ 

These formulae show not only that £ rapidly diverges as g — > 4+, but also that the 
leading behaviour of Eq. (2) accurate over a very wide range of q values. For 
example the correction term in Eq. (2) is still of the order of 1% for q as large as 75. 
The pure phase specific heats C Q and Ca are unknown, but their known difference 
vanishes when q — > 4 + as i"i These properties are of course in accordance with 
the point ( q = 4, j3 = fit ) being a second-order transition point, and lead one to 
speculate that Ca ~ C diverges as q — > 4+, possibly in such a way that the ratio 
C/£ is finite on both sides of q = 4. 

First we shall take advantage of the fact that, for the correlation length and 
latent heat, the "small q — 4 " region extends in practice up to large q values, 
and start from the opposite end. We compute the free energy of the model in 
the framework of a large q expansion, extrapolate down in q as far as we can, 
and analyze the resulting energy cumulants as functions of x in an intermediate 
q- value region. Nice regularities emerge, among which a smooth behaviour of C/x 
is ascertained. 

The large q expansion of the ordered free energy F Q (that for Fd in the dis- 
ordered phase follows from duality), was obtained through the Fortuin-Kasteleyn 
[6] representation of the Potts model partition function 

Z = - *) V (3) 

x 

where X is any configuration of bonds on a cubic lattice, / its number of bonds 
and n its number of connected components, or clusters of sites ( two sites bound to 
each other belong to the same cluster, an isolated site is a cluster). The completely 
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ordered configuration corresponds to n = 1 and / = 2V. So the partition function 
can be reorganized as an expansion in about this configuration: 



Z = q(eP-l) 2V £ Nl , n {V)^-J-)- l q n -i, (4) 

Z>0,n>0 ^ 



where Ni >n (V) is the number of configurations in a volume V with / removed bonds 
and n + 1 clusters. The enumeration of all the iVj jn (V) such that (/ — 2n) < M 
yields an expansion of Z Q to order M. Details will be given elsewhere [7]. To 
any given finite order M, a large enough volume V can be chosen to eliminate 
all boundary terms, so that all configurations retained correspond to disordered 
islands in a bulk ordered phase. We check that the sum in Eq. (4) exponentiates 
in V up to terms of order M + 1, defining a series for F Q truncated beyond order 
M. We have computed up to order M = 9, including terms up to N^g^oiV)- This 
series, whose first terms can be compared to existing low temperature series [8], 
provides us with similar series for the k th derivative with respect to ft, Fo h \ At 
P = f3 t the k = (free energy) and k = 1 (internal energy) series match the exact 
results [3] up to M = 9. The k = 2 and 3 cases give 

16 34 114 254 882 1944 6128 13550 
F ( ' = 1 1 1 1 1 1 1 (5) 

q ^ q 2 q 5/2 + q 3 ^ q 7/2 + q 4 q 9/2 ' 

/on 64 430 2654 12186 57018 224732 888024 3164682 , . 

_p(3) — Q) 

q ? 3/2 ? 2 ?5 /2 q 3 q 7/2 q 4 q 9/2 " ^ 

Let us make a few comments on these expansions. 

(i) F (2) gives the specific heat C Q = /3 t 2 F G (2) . The F Q (3) series gives - ((E - E ) 3 ), 
the first odd moment of the energy distribution in this phase. 

(ii) All terms in each series have the same sign, being of course positive while 

i?(3) < 

means that E > E Q is favoured with respect to E < E Q . 

(iii) The coefficients are fastly increasing with the order, the more so for larger 
values of k. This confirms the expectation that huge energy fluctuations 
are associated with the large correlation length, of order x in Eq. (2) when 
q decreases towards q = 4. This is expected to have a direct impact on 
the numerical analysis of these models, large values of the high cumulants 
invalidating the commonly used two gaussian formula [9]. 

To proceed with a quantitative analysis of F^ and F^ as functions of q, 
we conjecture for these quantities an essential singularity at q = 4 as £ has and 
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construct Pade approximants for the series of \nF^ instead of F^ k \ i.e. take as 
estimates of 

F^=e W [Pade(lnF^)} (7) 

The result of this construction for F^ and F^'/F^ 2 ' as functions of x for q = 
30, 20, 15, 10, 8, 7 and 6 is summarized in Fig. 1 as a log-log plot. The error bars 
are rough estimates of the uncertainties resulting from the higher order terms and 
have been obtained by varying the degrees of the numerator and denominator 
of the Pade approximant. The lines represent our prejudices F^ ~ Cst x and 
_p( 3 ) /_p( 2 ) ~ Cst x 3 / 2 (see below), where the constants are fixed by the q = 10 
values. It is clear that the general trend of both quantities is well reproduced over 
an astonishingly large range by such simple forms. 

As a by-product of this study we obtain analytical estimates for the ordered 
phase specific heat, which are compared to existing numerical data in Table 1. 



4 



Q 






Ref 


20 


5.362(3) 


5.2(2) 


[10] 


10 


18.06(4) 


12.7(3) 


[11] 






~ 18 


[11] 






10.7(1.0) 


[12] 


8 


37.5(4) 


23.(3.) 


[12] 


7 


71.3(1.0) 


47.5(2.5) 


[13] 






50.(10.) 


[11] 






44.4(2.2) 


[14] 



Table 1 

The ordered phase specific heat, our prediction compared to numerical estimates 

At q = 20 our prediction is in good agreement with the numerical estimate. 
In contrast it strongly disagrees at q < 10 with the value of C Q obtained from a 
finite size analysis of the maximum of the specific heat measured in the coexistence 
regime. However it agrees at q = 10 with the estimate obtained in [11] at (3 = /? t , 
in accordance with the rigorous statements of [1] . 

The large q expansion analysis supports the idea that not only does the cor- 
relation length in a pure phase diverge at (3 = /5 t , q — > 4+, but also that the 
associated fluctuations imply divergences of the energy cumulants. For example 
we obtain Fs ~ —1800 at q = 10 (see Fig. 1). Moreover the internal energy 
fluctuations behave in a way consistent with C/£ being finite at q = 4 + as it is 
known to be at q = 4_. We then propose the following ansatz: 

There exists a region of q > 4 where the free energies around f3 = f3 t reflect 
accurately the scaling properties associated with the second-order point lying 
at q = 4, (3 = In 3 and characterized by the corresponding critical indices a 
and v. 

Specifically, according to the known value 2/3 of a, we parametrize Fq 2 \(3) at 
q = 4 and for f3 — > (A)+ as 

F^((3) = A((3-(3 t )-V\ (8) 
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The higher derivatives are trivially deduced and in their expressions we replace 
(P — P t ) by ( Cst £~ 3 / 2 ) from Eq. (1) and, boldly continuing above q = 4 at 
P = fit, reexpress Fo P+2 ^ as a function of x via Eq. (2) to get 

F ( P+ 2) = Fo (h-2 )(a) = A ( _ )P r(2/3 + p) (Bg)1+3p/2i (g) 

where the constant S takes into account proportionality constants. Thus we get 
the following representation for the ordered phase free energy 

F ((3) = F(p t ) - E (P - p t ) + - Pt) n — - (10) 

n=2 

where we introduced the known linear term (F^ = —E ), and took F^ as given 
by Eq. (9) for n > 1. Note that all the odd cumulants are negative, as we found 
to be the case for F^ from the large q expansion. A similar expression holds for 
Fa, starting from Eq. (8) with (f3 — f3 t ) — > (fit — fi), and replacing E Q by Ed in 
Eq. (10). 

It is easy to sum the series Eq. (10). For later convenience, we introduce 
scaled temperature (v) : energy (e) and length (S) variables 

. = 0- . , S> = ^ (11) 

and end up with the following compact result 

FM = F(A) + ]p [~l - (eo + 1)« + |(1 + ^) 4/3 ] (12) 

This equation is the central result of this letter. We claim that, although we 
neglected all regular and less singular contributions to Eq. (8), Eq. (12) summarizes 
accurately, over a wide range of q > 4 values, all the properties of the model. 

Let us justify this statement by comparing the predictions of Eq. (12) for 
the energy distribution to data [11] taken at q = 10 with various L values. This 
distribution is obtained by inverse Laplace transform of the partition function 



/•/3 +«oo 

P V (E) = N t / dp qexp[VF (P)} + e W [VF d (P)] exp[VE(P - ft)] (13) 

Jfia—ioo 



where Fd follows from Eq. (12) by v — > — v and e — > ea (this is consistent with 
duality up to terms of order x~ 3 / 2 as compared to 1) and with N\ a factor ensuring 
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the probability normalization. Trading (/3 — /? t ) for v and E for e of Eq. (11), we 
get 

Pv(E) = N 2 / dve W [(-ne-e -l)v + -(l + vr/ 3 )]} (14) 

Jv —ioo & ^ 

Details on the computation of this integral will be given in [7] and here we limit 
ourselves to short remarks. 

i) As a function of e — e , Py depends on q and L only through the scaled volume 
(L/S) 2 . 

ii) Any vq > — 1 is suitable and the integral converges exponentially. 
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iii) At large (L/S), a saddle point method can be valuable. However (L/S) is 
not large in practice and there exists an energy value slightly above E Q where 
the saddle point value reaches v = — 1 (metastability point). 

iv) Actual computation requires numerical values of A and B. At q = 10 we get 
A = .193 and B = .386 from the large q expansion results (at other q values, 
slight changes have to be made according to Fig. 1 ). 

v) At q = 10 the length scale is S ~ 60, nearly 6 times the correlation length, so 
that for current values of L the ratio (L/S) is hardly of order 1 ! 

For the above reasons, we compute the integral Eq. (14) numerically. The 
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results are shown together with the data of [11] on Fig. 2 for L = 16, 20,24, and 
on Fig. 3 for L = 36, 44, 50. Remembering that the continuous curves are absolute 
predictions without any free parameter, it is quite striking to see how such a simple 
ansatz as Eq. (14) yields good results. Because (L/S) is not large, they are very 
different from what an asymptotic expansion would give. For example -Speak — E Q 
behaves effectively as ~ 1/L over a wide range of L values, whereas the asymptotic 
expectation (saddle point method in Eq. (13)) is 

L 2 (E - E peak ) = -Fo (3) /(2Fo (2) ) ~ 15, 100,350, 1000 at q = 20, 10, 8, 7. 

Marked discrepancies only appear at the external edges of the ordered (left) 
and disordered (right) peaks. In particular the theoretical curve levels out unduly 
at the bottom right of Fig. 2. This is a spurious effect : the tail of the ordered 
peak contributes more than the disordered phase. However it appears at a neg- 
ligibly small level at larger L's, being asymptotically of order exp[— c L 2 ]. As L 
is increased (Fig. 3), the agreement between predictions and data becomes better 
and better at nearly all values of E, but around the dip between the peaks. There 
indeed mixed phase contributions with percolating interfaces should finally win 
over the pure phase contributions. We consider the small departure of the theo- 
retical curve below the data points at L = 50 as an evidence for the emergence of 
mixed phase contributions. Since the dip region is often used for the determination 
of the interface tension a because strip configurations eventually yield an energy 
independent plateau in Py(E) with [15] 

2a- 1 In ^igdip) ^ (15) 

L Py (Speak)' 

the value and L dependence of the right hand side of Eq. (15) are interesting 
issues. Although in our construction this quantity diverges as L asymptotically, 
we unexpectedly find it roughly constant for L < 50, at a value around .11, not 
very far from the exact value l/£d = 0.95 [5]. This casts some doubt on attempts 
to determine a from Eq. (15) [13,14] unless the plateau in E is seen [16,17] 

Summarizing, we have shown, for the 2-d Potts model at q > 4, up to large 
q values, that the bulk properties are accurately described close to f3 t by free 
energies F Q {(3) and F d (f3) inferred by a simple ansatz from the known properties 
of the second-order point at q = 4. Definite expressions for F (f3) and F d (f3) 
were obtained by matching this ansatz to a calculation of their large q expansion, 
achieved at order (1/y/q) 9 . Our main result is expressed in Eq. (12) and illustrated 
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by the adequacy of the corresponding energy distribution to explain most of the 
features observed in numerical simulations at finite volume. 

In the first-order region, each phase knows little about the existence of the 
other ones, and rather feels the effectively close critical point q = 4, (3 = /? t , from 
which it inherits nice scaling properties. This is so as long as the linear size of 
the box is smaller than the length scale <S, proportional to the correlation length 
£, but many times larger ( S = 5 to 7£ in the range q = 7 to 20). Only at 
V > S 2 asymptotics takes place, so that interface tension effects can show up and 
be measured. 

Similar ideas could be applied to other situations. One such situation is the 
3-d Potts model at q > 3 where the first-order regime could also be influenced by 
the second-order point q S d situated between q = 2 (Ising) and q = 3. Then one 
might get some information on this critical point from numerical studies at q > 3. 
Another interesting case is QCD at finite temperature; although the transition is 
first-order, universality might be nevertheless invoked to relate its behaviour to 
that of the 3-d 3 states Potts model, both models sharing the same (universal) 
behaviour at a " close" point of parameter space. 

We are aware of the fact that our ansatz describes the pure phase free energies 
as analytic functions of j3 at /?t, which they probably are not, in the same way as 
in field driven first-order transitions the zero-field point is an essential singularity 
of the free energy [18]. This question deserves further discussion, but we believe 
that our ansatz, although not 'analytically correct', takes into account the most 
significant features of the model. 

It is a pleasure to thank R. Balian and A. Billoire for illuminating discussions. 
We also acknowledge useful conversations with B. Grossmann and T. Neuhaus. 
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